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HANKEL DETERMINANT SOLUTIONS TO SEVERAL DISCRETE 
INTEGRABLE SYSTEMS AND THE LAURENT PROPERTY 

XIANG-KE CHANG*? XING-BIAO HU? GUOCE XINt 

Abstract. Many discrete integrable systems exhibit the Laurent phenomenon. In this paper, we 
investigate three integrable systems: the Somos-4 recurrence, the Somos-5 recurrence and a system 
related to so-called Ai Q-system, whose general solutions are derived in terms of Hankel determinant. 
As a result, we directly confirm that they satisfy the Laurent property. Additionally, it is shown 
that the Somos-5 recurrence can be viewed as a specified Backlund transformation of the Somos-4 
recurrence. The related topics about Somos polynomials are also studied. 
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1. Introduction. Laurent phenomenon is a crucial property behind integrality 
shared by a class of combinatorial models while integrability is a key feature for a 
class of what we call integrable systems. An interesting observation is that many 
discrete integrable systems exhibit the Laurent phenomenon, and many mappings 
with the Laurent property are proved to be integrable. In the recent decade, a lot of 
related work were done via the study of a class of commutative algebras, called cluster 
algebras [niET], for its success in proving the Laurent phenomenon. For instance, 
see [m [121 [H 1221 ESI 1221 ESI ESI EH ES] • More recently, a generalization of cluster 
algebras called Laurent phenomenon algebras [39] was also introduced in order to 
study the Laurent phenomenon [H [JS] . 

One question is whether we can prove the Laurent property (or a stronger prop¬ 
erty) from other views. In this paper, we shall give one possible choice—the view of 
explicit determinant solution. That is, the Laurent property of some discrete inte¬ 
grable systems can be proved by their determinant solutions. Of course, how to obtain 
the desired determinant formulae is another challenging problem. This approach suc¬ 
ceeds for the following three integrable systems: 

The first two systems are the Somos-4 recurrence and the Somos-5 recurrence 
defined by 


SnSn-4 = aSn-lS„-3 + ^S^_2, ( 1 . 1 ) 

and 

SnSn-5 = dSn-lSn-4 + PSn-2Sn-3, ( 1 - 2 ) 

with nonzero parameters a,/3,Q;,/3 and arbitrary initial values, respectively. Here we 
assume Si never vanish. By using their theory of cluster algebras [B EH, Fomin 
and Zelevinsky [B E21 ES] proved that the Somos-4 and Somos-5 recurrences exhibit 
the Laurent property, that is, 5'„ are Laurent polynomials in its initial values with 
coefficients in Z[a,,5] and Z[q:,/ 3], respectively. Hone constructed explicit solutions 
to the Somos-4 and Somos-5 recurrences in terms of the Weierstrass sigma function 
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in [3D] and [33], respectively. He also indicated that the Somos-4 recurrence can be 
thought of as an integrable symplectic mapping [3115] by making a change of variables 
and Somos-5 be a integrable Poisson mapping. Unfortunately, it is not obvious to see 
the Laurent phenomenon from their Weierstrass sigma function formulae. 

Our result for the Somos-4 recurrence with initial values 1,1, x^y is as follows. 

Theorem 1.1. If we let S-i = l,S'o = IjS*! = x,S 2 = y,S 3 = ay + fix"^ 
by shifting the indices, then the Somos-4 recurrence (HH) has the following explicit 
Hankel determinant representation: 





PkPm- 2 -k{m > 2) with 


where Pm 


initial values po = x,pi = —yfa. 

The Laurent phenomenon clearly holds for initial values 1,1, cc, y by the theorem: 
i) the determinant formula asserts that Sn is a Laurent polynomial in ^/a, x, y and 
polynomial in /3; ii) the recursion (inD shows that Sn is rational in a and that Sn 
exists and not equal to 0 when a = 0, which can be checked by induction. To see 
that the Laurent phenomenon holds for arbitrary initial values, we observe that the 
recursion holds for Sn if and only if it holds for 5'„/c for constant c^Q, if and 
only if it holds for Snf^ for constant t ^ 0. This makes it sufficient to consider the 
sequence {f^(f 7 )"}n>o, which starts with 1,1. Note that the reduction to the case 
1,1, x,y is a simple but important step, since the Somos-4 recurrence with arbitrary 
initial values does not seem to have the desired Hankel determinant representation. 

Similarly, it is sufficient to give the following result for the Somos-5 recurrence. 

Theorem 1.2. If we let S -2 = 1, S-i = 1, S'o = x,Si = y, S 2 = z, S 3 = az-\-j3xy, 
S 4 = cP'xz + afix'^y Pyz by shifting the indices, then the Somos-5 recurrence (lEa 
has the following explicit Hankel determinant representation: 


^271 — ^ — 

5*2n+l — y ci6t )o< 2 ,j<n —1; 
where Pm is recursively given by pq = 


(1.3) 

(1.4) 


—ax'^z -\- ax'^y'^ -\- jdx^y -\- zy"^ — z^ 



and qm is recursively given by go = = —P, 


ax^z ax^y^ — z^ — zy^ -|- Px^y 


Qm — 


Qm—1 


xyz 



m—2 


The third system is called the extended Ai Q-system: 

SnSn—2 = Sn-i + P, 


(1.5) 
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with nonzero parameters /3 and arbitrary initial values. Here we assume Si never 
vanish. The special case when (3=1 and all the initial values are Is reduces to the Ai 
Q-system [HIIlIIl]. From the view of integrable systems, (11.51) is C-integrable[5]. 
Note that the Laurent property of m was shown in [niiioiiM]. 

By dividing (11.51) by S'q, we may assume the initial value is S'o = 1, S'! = x (/3 
changes but do not affect the Laurent property). The determinant solution is as 
follows, from which the Laurent property is obviously satisfied. 

Theorem 1.3. If we let Sq = l,5'i = x, then the extended Ai Q-system (Hm 
has the following explicit Hankel determinant representation: 

Sn — det (pi+j )o<zy <n—1; 

where Pm = Pm-i + PkPm-i-ki'm > 2) with initial values po = x,pi = 

1 . 

All three theorems can be proved by employing Sulanke and Xin’s method of 
quadratic transformation for Hankel determinants [43l [48] or classic determinant tech¬ 
nique. In particular, the special case x = y = a = /3 = loi Theorem 11.11 was obtained 
by the third named author in [48] . It is also confirmed by using the classic determinant 
technique in [^, which is used to solve a conjecture that a certain determinant satisfies 
a specified Somos-4 recurrence US]- In Section 3, we shall give a detailed proof for 
the Somos-4 case by using Sulanke and Xin’s method of quadratic transformation, 
and prove the result for the extended Ai Q-system by applying the classic determi¬ 
nant technique. As for the Somos-5 recurrence, we can solve it via the Backlund 
transformation |291 140] of the Somos-4 recurrence, which will be introduced in Section 
2. As a result, the Laurent property for the three models are confirmed from their 
determinant solutions. In Section 4, we prove the Somos-4 polynomials result in [41] 
in our own way and give more Somos polynomial sequences. 


2. Preliminaries. 

2.1. Sulanke and Xin’s quadratic transformation for Hankel determi¬ 
nants. There are many classical tools for evaluating Hankel determinants, such as 
orthogonal polynomials approach (cf. e.g.[38l[Tni|4^), Gessel-Viennot-Lindstrom the¬ 
orem (cf. ep.g. O [27l |43]),the J -fractions (cf. e.g. [38l|47]) and the S-fractions (cf. 
e.g.[37]). Recently, Sulanke and Xin proposed a method of quadratic transformation 
for Hankel determinants [43] developed from the continued fraction method of Ges- 
sel and Xin |5S]. In [35], Xin applied the special quadratic transformation to solve 
Somos’ conjecture about the Hankel determinant solution to the Somos-4 recurrence 
with initial values Si = l,i = 0,1, 2, 3 and a = l,/3 = 1. We restate the iterative 
transformation as follows. 

Lemma 2.1. Given the initial values ao,bo,co,do,eo, fo, let the generating func¬ 
tion Qo(x) be the unique power series solution of 


Qo{x) 


ao + box 

1 -I- cox dox"^ + a;2(eo -I- fox)Qo{x) 
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//a„+i,6„+i,c„+i,(i„+i,e„+i,/„+i are recursively defined by 


^n+1 — 




bn+l — -- 
^n+l — 
dn-\-l = — 


^fn “t“ CnCL^dji CL^Cj^bji -\- ‘2ciYiCnb^ Cl^bfidn b^ 


drf dji ^OinbnCn ^b^ 


^n+1 I5 

P _ 

/n +1 5 

dji 

where we suppose that ^ 0, n = 0,1, 2, • • • ^ then 

det{Hn{Qo)) = Ooa"”^ • • • a„_i. 


Here we remark that Hn{R) denotes the Hankel determinant det(ri+j)o<z,j<n-i for 
any power series R{x) = 

The basic tool in our approach is the following surprising result, which solves the 
above recursion system. 

Lemma 2.2 (Theorem 2 in [48]). Suppose Cn = c,en = —1, then the recursion 
system in Lemma\2f^ satisfies 


0n+20n+l + On+lOn — 2aoai + Oq (2 /i + c) (/o + C + ff) 


ap^ ifo + c + fif 

^n+1 


2.2. Backlund transformation. A Backlund transformation is a transforma¬ 
tion between a solution m of a given differential or difference equation, 


£i{u) = 0, 

and another solution v of another differential or difference equation, 


Mv) = 0, 


where £2 may be the same as, or different from, Ci. It is an important tool for finding 
new solutions in soliton theory and integrable systems. See [29l [40] for more details. 
Here we shall give an example, which is concerned about Backlund transformation 
in bilinear form between two solutions of the Somos-4 recurrence. It can be used for 
solving the Somos-5 recurrence. 

Firstly, we rewrite the Somos-4 recurrence as a bilinear form 

- f)S„ ■ S„ = 0, (2.2) 


where the Hirota’s bilinear operators [55] are defined as follows: 

or ■ 0= tw. 01..=.,,.=,, 

d d 

e^^'^Qn ■ bn = exp [H-^ - -^)]anbn'\n'=n = Qn+sbu-S- 
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Then, we have the following lemma. 

Lemma 2.3. The Somos-4 recurrence has the following bilinear Bdcklund trans¬ 
formation: 


(e - Ae^" - fi)fri ■ gn = 0, (2.3) 

- ae^^^^"-)fn • = 0, (2.4) 

where /t, A, ry are arbitrary constants. 

Proof Let /„ be a solution of the Somos-4 recurrence. What we need to prove is 
that gn satisfying (lOD and (12.41) is another solution of the Somos-4 recurrence, i.e., 

P = - /3)g„ • 5n = 0, 

Then, using the bilinear operator identities 

{e^’^'fn ■ fn)gl - fl(.e^^’'gn ' 5n) = 2sinh(D„)/„ • g^) • fn • 5n), 

(e^"/n • fn)gl - fl{e^-gn ' gn) = 2 sinh(l/2i4„) /n • 9n) ■ {e-^^^^-fn ■ gn), 

sinh(D„)(e^"/„ • gn) ■ (fn ■ gn) = sinh(l/2D„)(e^/^^'‘/„ • gn) ■ • gn), 

smh{5Dn)fn ■ fn=0, 

we have 

-Pf^ = - P)fn • M - flie^^- - ae^- - /?)g„ • 

= [i.e^’^-fn • fn)gl - fl{e^^-gn ' 9n)\ - a[(e^"/„ • fn)gl - fl{e^-9n • gn)\ 

= 2sinh(D„)(e^'‘/„ • gn) ■ {e~°^fn ■ gn) 

-2asinh(l/27A„)(ei/2^"/„ • gn) ■ • g„) 

= 2/rsinh(L)„)(e-°’‘/„ • gn) ■ [fn ■ gn) 

-2yisinh(l/2D„)(e^/^^"/„ • gn) ■ fn ■ gn) 

= 0 . 


Thus the proof is completed. □ 

3. Determinant solution and Laurent property. 

3.1. On the Somos-4 recursion system. The Hankel determinant formula in 
Theorem o is discovered by employing Sulanke and Xin’s quadratic transformation 
for Hankel determinants. The proof is given below. 

Proof of Theorem M.R Let P{t) = generating function of the 

sequence {pn}^=o^ need to prove S'„ = det(iL„(P)). Additionally, it is easy 

to see that P{t) satisfies 


m 


_«o + bpt _ 

1 -I- Cot -I- dof^ + t2(eo -I- fot)P{t) 


(3.1) 
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where 


Q-0 — ^5 


bn = — 


Co = 


dn — — 


+ ax^ 

y/ay 

/3x‘^ — y^ + ax^ — ay 
yfaxy 
j3 + ax — xy 


eo = -1, 

/o = 0. 


By Lemma l2.11 the recursion m we need to prove is transformed to 


^n^n—l^n—2 O: “h f (^n — 1- 

r.2 


(3.2) 


Lemma [2.11 also gives oi = y/x^ and fi = ' Thus applying Lemma 

2.21 gives 


Px'^ + y'^ + ax^ + ay ^ 

C'n+2 — ( j/C'n+1 ^/C-n+l C-n* 

xy 


(3.3) 


By substituting (13.31) with n replaced by n — 2 into (13.21) , we need to show that 


T[n) := ( 


Px"^ + + ax^ + ay 


xy 


)a„_ia„_2 - aa„_2 - - aa„_i - P 


0 (3.4) 


holds for n > 2. 

We prove this by induction. It is easy to confirm T{2) = 0. Assume T{n — 1) = 0. 
Replacing an-i by (1531) with n replaced by n — 3 in T{n), we have 


Tin) = ( 


Px^ + + ax^ + ay 

xy 


)a„_2a„-3-Q!an-3-a^-2an-3“««"-2-/3 = T{n-1) = 0. 


This completes the proof. □ 

3.2. On the Somos-5 recursion system. In finding the determinant solution 
of the Somos-5 recursion, the key observation is that the even and odd terms of the 
Somos-5 sequence are both Somos-4 sequence. This observation has been made by 
Hone (Proposition 2.8 in [33]), through making the connection with a second order 
nonlinear mapping with a first integral. Here we give an alternative explanation from 
the view of Backlund transformation, which seems simpler and more intuitive. Now 
we prove Theorem 11.21 by using Lemma 12.31 and Theorem 11.11 

Proof of Theorem II.^l If we let /„ = S 2 n and y„ = 5271 - 1-15 then we get two 
coupled equations 


9n+lfn—l d^fn+l9n — l Pdnfn — 0 ; (^■^) 

fn+2gn-l ~ agn+lfn ~ Pfn+ign = 0. (3-6) 

Recall that, as is indicated in Lemma [2.31 the Somos-4 recurrence dm has the 
following Backlund transformation 
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9'n-\-lfn—l -^/n+l^n—1 99nfn — Oj 

9'fn-\-29n—l '99n-\-lfn (^fn-\-l9n ~ Oj 


(3.7) 

(3.8) 


where /i, A, ij are arbitrary constants. When fj, = p, X = a, ij = ap, a = /3^, it is easy to 
check that (13.71) . (13.81) are the same as (13.5L(13.61) . This means that subsequences /„ = 
S 2 n and gn = S 2 n+i of the Somos-5 recurrence dm) may satisfy the same Somos-4 re¬ 
currence dm with a = /32 and /3 = 

xy z 

where /3 is determined by the first five values of /„. In fact, by use of the meaning 
of Backlund transformation, this can be confirmed by induction. Thus, it is easy to 
construct the general solution to the Somos-5 recurrence dm- And one can complete 
the proof after some substitutions and calculations. □ 

It is natural to ask if we can give a more direct proof to the Somos-4 recursion of 
S 2 n (or S 2 n-i)- The answer is positive. We sketch the idea for S 2 n in the end of this 
subsection. It is easy to perform the detailed proof by Maple. 

We can eliminate the odd Sn to obtain the recursion 


S2n-&S2n — >^2 n-4*2 n-2 


S 2 n -2 


/3 (S'2n-6<S'2n-2 + d S'2n-4^) 


^ 2 n— 8 ^ 2 n —2 


P (S'2n-6'S'2n-2 + d 5'2n-4^) 
using which we can prove 

^2n^2n—8 — OiS2n-2S2n-& + /3'S'2n-4 


^2 n—6*^2 n—4 


(3.9) 


by induction. 

Eliminating S 2 n in the above equation by using ()3.9I) . we are left to show the 
nullity of a relation, saying target{S 2 n- 2 , <S' 2 n- 4 , S 2 n- 6 , S 2 n- 8 )- Eliminating S 2 n -2 in 
target{S 2 n- 2 , S 2 n- 4 , >S' 2 n- 6 j <S' 2 n- 8 ) by using (13.9p with n replaced by n — 1, we obtain 


tCLrg€t(^S2n—2^ 5*2n—4i S2n—Qi ^271—s') — 


‘^2n—4'^2n—10 
‘^2n-6'^2n-8 


t(lTg€t(^S2n—4:: *5*2n—6; *5*2n—8; »5*2n —lo) 


Then the proof is completed after checking target{Se, S 4 , 82 , 80 ) = 0. 

3.3. On the extended Ai Q-system. Let us first prove Theorem [LS] by classic 
determinant technique. This approach gives no hint how we discovered the Hankel 
determinant formula, but this proof seems to be shorter than using the Sulanke and 
Xin’s quadratic transformation method. 

Proof of Theorem, \1.S[ Let 77® denote det{pi+j+i)o<i,j<n-i and we shall use the 
conventions 

_ I 

m ’ (3-10) 

H):> = 0 forn < 0. 

Firstly, we assert that there hold 

77« = l, n>0; (3.11) 

77 ( 0 ) = + pH^22, n>l. (3.12) 

Employing the Jacobi determinant identity HIS], we get 

77^77^2^ = - iH^nli)^ ri>2. (3.13) 
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By use of (13.lip and (I3.12L replace in p.l3p . then we can 

obtain 


= + n>l. 

Noting that S'„ = Hn \ thus it suffices to confirm the formulae (13.111) and (13.121) . 
It is noted that Pm also satisfy the following recurrence: 


+ 1 


Pm = - Pm-1 + PkPm-l-k, TTl > 2. 

nr. 


(3.14) 


fe=i 


We will prove (13.111) and (I3.12p by employing row and column operations for the 
determinants and using this recurrence. 

Let’s consider firstly. 

Step 1 : Subtract the i-th column multiplied by Pn-i-i from the n-th column for 
i = 1,2,'' - ,n — 1 and subtract the (n — l)-th column multiplied by from the 
n-th column. Next, applying a similar procedure to the (n— l)-th,(n — 2)-th,- • • , 2-nd 
columns, by using the recursion relation (I3.14p . we have 



Pi 

0 

0 

P2 

PiPi 

2 

PlPn-1 

2 

P3 


... '^PiPn+l- 


2=1 

2=1 


n — 1 

n—1 

Pn 

E PiPri-^ 
2=1 

... ^ PtP2n-2 

2 = 1 

PlPl 

2 

PlPn-1 

2 

E 

2=1 

PiP3-i ■■■ 

^ ^ PiPn+l—i 

2=1 


n —1 

X! P^Pn-i 


n—1 

y^PiP2n-2-i 

i=l 


Step 2: Perform row operations for the above determinant. For fixed k = 
2, 3, • • • , n — 1, we subtract the i-th row multiplied by pk+i-i/pi for f = 1, • • • , fc — 1 
from the fc—th row. Then it follows that (13.111) holds. 

Now let’s turn to the proof of (13.121) . Obviously, the result holds for n = 1. In 
the following we consider the case of n > I. 

Step 1: Subtract the i-th column multiplied by Pn-i-i from the n-th column for 
i = 2,3,-'' ,n — I and subtract the (n — l)-th column multiplied by from the 
n-th column. Next, applying a similar procedure to the (n— l)-th,(n — 2)-th,- • • , 2-nd 
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columns, by using the recursion relation (I3.14|) . we have 




PQ 

-P 

0 

0 

Pi 

PlPO 

2 

PlPl 

2 

PlPn-2 

2 

P2 

yypiP2-z 

yypzP3-i ■ 

y^PzPn-i 


2=1 

2=1 

2=1 


n—1 

n—1 

n—1 

Pn-1 

^ ^ PiPn — l — i 
2=1 

yy PiPn-z ■ 

2 = 1 

y^^PiP2n-3-i 
2 = 1 


= X 


PlPO PlPl 

2 2 

^PtP2-t '^PiPS-z 


n—1 n—1 

^ ^ PiPn—l—i ^ ^ PiPn—i 
2=1 2=1 


PlPn-2 

2 

'^PiPn-i 

2=1 

n—1 

PiP2n-3-i 

i=l 


Pi 


+/3 


P2 


PlPl 

2 

i=l 


PlPn-2 

2 

y^PiPn-z 

i=l 


n—1 

Prz — l yyPiPn—i 
z=l 


n—1 

yyPiP2n-3-i 

i=l 


Step 2: Perform row operations for the above two determinants. For fixed k = 
2, 3, • • • , n — 1, we subtract the i-th row multiplied by pk+i-i/pi for t = 1, • • • , fc — 1 
from the fc—th row. Then it follows that (13.121) holds. 

Therefore, we complete the proof. □ 

Remark: If we let /„ = then /„ satisfy fnfn -2 = fn-i + with 

/o = 1 and fi = Xy/j. And we can see that {/«} gives the Fibonacci sequence when 
7 = —1, /3 = 1 and x = It means that every Fibonacci number can be 

expressed as a Hankel determinant. 

Moreover, it is noted that the extended Ai Q-system satisfy a three term re¬ 
currence [34) . In the following we use the relations of Hankel determinants to show 
it. 

Theorem 3.1. The extended Ai Q-system with Sq = 1, Si = x satisfy the 
following three term recurrence 

Sn = - — - ^^Sn-l-Sn-2- (3.15) 


Proof. We shall firstly prove the formula 

= xH^^y - 


(3.16) 
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where the notations are the same as that in the previous proof. It is easy to see that 
Hn'^ = + Hn, where 


0 

Pi 

P2 

■ • Pn-1 

Pi 

P2 

P3 

Pn 

P2 

P3 

P4 

Pn+1 

Pn-1 

Pn 

Pn+1 ■ ' 

P2n-1 


( 2 ) 

Based on this observation, it suffices to prove that Hn = — 

Step 1: Subtract the i-th column multiplied by Pn-i-i from the n-th column for 
i = - ,n — 1 and subtract the (n — l)-th column multiplied by from the 

n-th column. Next, applying a similar procedure to the (n— l)-th,(n — 2)-th,- • • , 3-rd 
columns, by using the recursion relation (I3.14[) . we have 


Hn 


0 

Pi 

0 

Pi 

P2 

0 

P2 

P3 

P 1 P 2 

2 

P3 

Pi 

'^PtPi-t 






n-2 

Pn-1 

Pn 

X! P^Pn-i 


i=l 


PlP2 

2 

'^PiPA-i 

n — 2 

X! PiPn-i 

i=l 


0 
0 

PlPn-1 

2 

^ ^ PiPn+l—i 

n—2 

■■■ y^PiP2n-3-i 

i=l 

PlPn-1 

2 

^ ^ PiPn+l—i 
i=l 

n-2 

PiP2n-3-i 

i=l 


Step 2: Perform row operations for the above determinant. For fixed k = 

2, 3, • • • , n — 2, we subtract the i-th row multiplied by pk+i-i/pi for z = 1, • • • ,k — 1 

~ (2) 

from the fc—th row. Then it follows that Hn = — 

Next, noting that Sn = Hn\ the three term recurrence can be obtained by 
combining p.l2l) and (13.161) . □ 

Remark: When /3 = — 1, {5'„} satisfy the three term recurrence 

Sn — 2x5*71—1 Sn—2- 


with initial data 5o = 1 and 5i = x. This just meets the Chebyshev polynomials. 

It is also noted that every sequence produced by the extended Ai Q-system is a 
Somos-4 sequence [44]. More precisely, if { 5 , 7 } satisfy (II.5L then { 5 , 7 } also satisfy 


5t7-|-2577_2 - 


(x^ 


■^)% o ,.7 (^^ + i+/3)yo2 
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This can be easily confirmed by employing the three term recurrence. 

We conclude this subsection by sketch the idea of the alternative proof of Theorem 
11.31 It turns out to be a long journey to prove by Sulanke and Xin’s method and new 
transformations are needed. The proof is too lengthy to be typed here but is easy to 
perform using Maple. 

The generating function F{t) = pq + pit + ■ ■ ■ is uniquely determined by 

This simple form has to be transformed so that Lemma 12.21 applies. Indeed, by using 
Proposition 4.2 in [43], we obtain 


G{t) = 


b+x^ 


_O-f 

I' 


1 - + 2 + t 2 (_i + At) G{t)' 


with the relation Hn{F{t)) = a;"iL„_i(G(t)). The recursion (II.5p we need to show 
becomes 

r„r„_2 - = o, (3.17) 


where r„ = i7„(G(<)). 

Now we prove by induction the recursion (13.171) with n replaced by n+3. Eliminate 
Tn+3 by the following recursion 

Tn+sPn r„+2Pn-l + 2 b + 2 bx‘^ + 4 + 1 + (6 + 1 + x") 

Tn+ 2 Tn+l Tn+lTn X^ x'^Tn+lTn 

which is obtained by application of Lemma 12.21 

Successively eliminate T„+ 2 )T'n+i by recursion (13.171) with n replaced by n + 2, 
n + 1. We are left to show the nullity of a polynomial, say targetiTn, T^-i), of degree 
8 in Tn and degree 4 in Tn-i- Now eliminate T„ as before, we obtain a new polynomial 
which has T„_ 2 ) as a factor. The proof is then completed after checking 

targetiTi.To) = 0 . 

4. Some Somos-polynomials. Our starting point is to try to give a proof of 
the following statement of Somos in [41] by using the above result. 

Theorem 4.1 (Somos-4 Polynomials). The sequences produced by 

^ _xyzSn-iSn-3FxyzSl_2 

*^n—4 

with So = X, Si = 1, S 2 = ^, So = y are all polynomials in x,y, z. 

We obtain the following corollary as a consequence of Theorem [O] by making the 
appropriate substitutions and calculations. 

Corollary 4.2. If we let S -2 = x,S-i = Ij^o = l,5'i = y,S 2 = y^z + yz,So = 
xy^z'^ + xy^z + xy'^z'^ by shifting the indices, then any Somos-4 polynomial can be 
expressed as 


Sn — det (pi+j )o<z j <n—1; 
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where Pm = + (a; - y)Pm -2 + Yl'k=o PkPm- 2 -k{m > 2) with initial 

values po = y,pi = —y/xyz. 

Our first hope was to obtain nicer determinant with polynomial entries by some 
transformations, but we find a simple proof by just using the Laurent property. 

Proof. [Proof of Theorem l4.1| By the Laurent property, Sn are Laurent polynomial 
with denominator factors x and y. To see that x and y disappear in the denominator, 
we extend the recursion m and obtain 

S'-! = xz{x + 1 ), So = x,Si = 1 , S 2 = 1, S '3 = y, S 4 = yz{y + 1 ). 


Using the Laurent property for we see that the factor y is not in the 

denominator; Using the Laurent property for {5'„}ji>i, we see that the factor x is not 
in the denominator. This completes the proof. 0 

The idea of the proof suggests the existence of a stronger version of the Laurent 
property. Indeed such a version was proposed by Hone and Swart [35l Theorem 3.1]. 

Theorem 4.3. For the Somos-4 recurrence HU, every Sn belongs to Z[a, fi, a^ + 
I3T, S^^, Si, S 2 , S 3 ], where T is a quantity independent of n defined by 


T = 


SnSn+3 

^n+lSn+2 


a{ 


02 

SnSn+2 


g 2 

‘^n+2 


5 * 71+1 5 * 71+3 


-)+P- 


^n+lSn +2 
Sn Sn+3 


(4.2) 


With this theorem, we can obtain more general polynomial sequences. 

Theorem 4.4. The sequences produced by 

^ xyzwSn-iSn-3 + xyzwSf _2 

On — ^ 14.Oj 

with the initial values for the following cases all yield polynomials in x,y,z,w: (1): 
l,x,w,y; (2): x,l,w,y; (3): x,w,l,y; (4): x,w,y,l. 

Proof Noting that /3 = xyzw and T is independent of n, from the expression of 
T, we conclude that /3T must be polynomials for all the four cases. And we always 
shift the indices to make 5o = 1 without changing the value of T. Thus, the result 
follows by Theorem 14.31 □ 

Remark: By setting lu = 1, the special case (2) of the above theorem reduces to 
Theorem HU In this case as in Corollary 14.21 we have a = fd = xyzw, 5o = 1 and 
T = yz + xy + xz + z. 

Combining the idea of our proof of Theorem 14.11 and using Theorem 14.31 we can 
obtain more general result. 

Theorem 4.5. Let {xi}^i be a finite set of indeterminates and T be defined by 
(HU- Suppose { 577 } satisfy the recursion system (II.ip and a, j3, o^+fiT are polynomi¬ 
als in all Xi- If there exists a nonnegative integer r such that gcd(5o, 5i,..., Sr) = I 
and Sk are polynomials in all Xi for 0 < k < r + 3, then Sn for n > 0 are all 
polynomials in all Xi. 

Proof. We prove it by contradiction. 

Suppose that D is an irreducible polynomial appearing as a factor in the de¬ 
nominator of certain Sm- Apply Theorem 14.31 to the Somos-4 sequences {5,i}7i>^ for 
^ = 0,..., r. We see that D must be a factor of Si for each £. Thus D divides 
gcd(5o,..., Sr) = 1. A contradiction. □ 

Hone and Swart also gave a strong Laurent property for Somos-5 [351 Theorem 

3.7]. 
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Theorem 4.6. For the Somos-5 recurrence (11.21) . any {S'„} G Z[q;,/3,/3 + aT, 
, S 2 , S^, S 4 ], where T is a quantity independent of n defined by 


T = 


^n+lSn+2 


^n+l^n+4 , ~ / ^n+l^n+2 

+ a{ 


Sn+2Sn+3 


Sn Sn+3 1 Sn4-4: 


^n+2 ^n+3 ^ ^ ^n+2 


SnSn+4: 


(4.4) 


Consider the recurrence 

WXyzSn-lSn-4 + WXyzSn-2Sn-3 , , r\ 

Sn = -^-. (4.5) 

^n—5 

Similar to the proof of Theorem 14.41 we can also confirm the following result for 
Somos-5 polynomials. 

Theorem 4.7. The sequences produced by (ITSI) with the initial values for the 
following cases all yield polynomials inx,y,z,w: (1): 1, x, w, 1; J// (2)' x,l,w,y,l; 
(3): l,x,w,y,l; (4): l,x,l,w,y; (5): x,w,l,y,l; (6): x,l,w,l,y; (7): l,l,x,w,y; 
(8): X, (9): x,w,l,l,y; (10): x,w,y,1,1. 

Proof Here we only give a detailed proof for case (1) because the proofs for the 
other cases are similar. 

Noting that a = xyzw and T is independent of n, from the expression of T, we 
conclude that aT must be polynomials for all the four cases. Thus, from Theorem l4.6l 
and S'o = 1, 5”! = x, it is obvious that Sn is a Laurent polynomial in x and polynomial 
in y,z,w. If we shift the indices to make Sq = w, Si = 1 without changing the value 
of T, then Sn is also a Laurent polynomial in w and polynomial in x, y, z. Therefore, 
Sn must be a polynomial in x, y, z, w and the proof is completed. □ 

Similar to Theorem 14.51 using Theorem 14.61 we also have the following general 
result for Somos-5 polynomials. Here the details of proof are omitted. 

Theorem 4.8. Let {xi}fLi be a finite set of indeterminates and T be defined 
by (14.41) . Suppose {5'„} satisfy the recursion system (11.21) and a,j5, fi -\- aT are 
polynomials in all Xi. If there exists a positive integer r such that gcd{SoSi, S 1 S 2 , 

..., SrSr+i) = 1 and Sk are polynomials in all Xi for 0 < k < r + 4, then Sn for n >0 
are all polynomials in all Xi. 

5. Conclusion and discussions. We have derived determinant solutions to 
three discrete integrable systems and confirm that the three discrete systems exhibit 
the Laurent property. The simplest case is the extend Ai Q-system. In particular, we 
prove that every Fibonacci number can be expressed as a Hankel determinant. The 
general solution to the Somos-4 recurrence is given in terms of Hankel determinants. 
Observing that the Somos-5 recurrence can be viewed as a specihed Backlund trans¬ 
formation of the Somos-4 recurrence, we also derive the Hankel determinants solution 
to the Somos-5 recurrence. 

As for the Laurent property of the Somos-4,5 recurrences, it can also be proved 
by using the combinatorial models as reduction cases of the octahedron recurrence 
[42] (also appearing as T-system in [lUdl]). The Ai Q-system, as a special case of 
the Q-system, was also done in M- But the combinatorial method seems somewhat 
awkward to work with in practice. Our determinant formula consists of elements 
with convolution recurrence, and it is a closed form, which appears clearer and more 
intuitive. (Although the determinant solution to T system with arbitrary initial values 
was given in [HI [17], it is not a easy thing to construct the solutions to Somos-4,5 
recurrences as reduction cases.) It should be noted that the combinatorial models 
imply the positive Laurent property (that is, the nonzero coefficients of the Laurent 
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polynomials are all positive), while our determinant solutions do not seem to imply 
this property. The positive Laurent property is also an interesting topic. See [Hill 
etc. 

Finally, we shall describe how we discovered the Hankel determinants formulae. 
We explain it by taking Theorem 11.11 as an example. Actually, by applying Lemma 
[2T] and 12.21 we can solve the Somos-4 recurrence for Hankel determinant solution, 
and the expression of the solution is not unique. This is equivalent to solving for 
ao,bo,co,dQ,fo (noting that e„ = -1). 

Let us sketch the idea as follows. As before, shift the indices so that S-i = Sq = 
I,!?! = x,S 2 = y, 3^ = ay + . It is easy to see that ao = x and oi = y/x^. Similar 

to the proof to Theorem ll.il in order to make the Somos-4 recurrence hold, we need 
to equivalently make another formula {T(n) = 0) hold, which can be obtained by 
applying Lemma 12.21 Then we need the induction step T{n) — T{n — 1) = 0 and 
T(2) = 0, from which there needs to make sure that 

r aofo + aoco - k^/a , 

fi = -, fc = ±1 (5.1) 

Oo 

and 


alaj -f 2 agOi/o/i -f 3 almcofi + 2 alaif^ + alaicfo + alaiCo'^ - aj^/o 
-2Oq/oCo - 2Og/o/i - agCo^ - 2 alcofi - - p - aai = 0, (5.2) 

respectively. 

Additionally, from the recurrence relations, we have 
_ aleo + aldo - aoboco + bl 

Qj\ — ^ 

^0 

Noting that fi = —^ and oq, eg, ai is known, we can express 5o, cq, do in terms of 
a free number /o by using (ED, (ED and (ED- Thus, we have the following result. 

Theorem 5.1. If we let S-i = IjS'o = l,S'i = x,S 2 = y,S 3 = ay + fix'^ 
by shifting the indices, then the Somos-4 recurrence ED has the following explicit 
Hankel determinant representation: 


3n — det(pj+j )o<z j'<n—1 ■ 


Here sequence {pn} has the generating function 


m 


_Qq + bpt _ 

1 -f Cot dot"^ -f t^(eo + fot)P{f) ’ 


ao = X , 


bo = - 


Co = -- 


yfoV^x -\-j5x^ — y'^ — ay-\- x^a kay 
yy/a 

j5x^ — y'^ — ay x^a 2 yfoyfax 
xyy/a 


do — — 


•^—^[—fix^fo -f fix^ky/a ayfox — x^afo + x^a^^^k y^fox 


x'^yy/a 


—a^/'^yk — yfly/ax^ -f k^a^^'^y — x^yy/a y^y/a — y^ky/a], 
Co = -1 


with 
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and fc = ±1 and arbitrary number /q. 

Remark: The special case of A: = 1 and /o = 0 reduces to Theorem 11.11 which 
seems to be the most concise version. 

Additionally, it is noted that one can find artificial interpretations in terms of 
weighted sums of some combinatorial objects for the general Somos-4,5 and the ex¬ 
tended Ai Q-system by employing the well-known Gessel-Viennot-Lindstrom theorem 
(cf. e.g. [5l[2Illl3]) for determinants. 
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